Kink-antikink vortex transfer in periodic-plus-random pinning potential: 
Theoretical analysis and numerical experiments 
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The influence of random pinning on the vortex dynamics in a periodic square potential under 
an external drive is investigated. Using theoretical approach and numerical experiments, we found 
several dynamical phases of vortex motion that are different from the ones for a regular pinning 
potential. Vortex transfer is controlled by kinks and antikinks, which either preexist in the system 
or appear spontaneously in pairs and then propagate in groups. When kinks and antikinks collide, 
they annihilate. 
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The behavior of an elastic media under the competitive 
action of a regular potential and disorder is a common 
problem in various fields of modern physics. Examples 
of such media are vortex lattices in superconductors^*^ 
or in Bose-Einstein condensates of ultracold atoms^, in- 
teracting colloids on periodic substrates&i, charge and 
spin density waves in metals^, polarization density waves 
in ferroelectrics and many others. Regular pinning po- 
tential can be either of artificial origin, as in nanostruc- 
tured superconductors and in Bose-Einstein condensates 
with optical lattices, or it can be imposed by the crystal 
structure of the material. In superconductors, pinning ef- 
ficiency determines the value of the critical current which 
is of great practical importance. Theoretical description 
of such systems is a quite complicated problem, which in 
one-dimensional case can be reduced to the well-known 
Frenkel-Kontorova model^. 

Recently, two-dimensional lattice of repelling particles 
in the presence of square pinning potential and disorder 
was investigated theoretically and numerically 10 . In the 
case of weak disorder, a pinned vortex lattice is disturbed 
by specific defects consisting of elastic strings of depinned 
vortices. It was found that these strings are able to in- 
tersect and form branched fractal-like clusters, which can 
percolate through the system. The aim of the present let- 
ter is to investigate the vortex dynamics in such a system 
and look for correlations between the static defects and 
the dynamical phases under an external drive. 

Model. — We model a three-dimensional (3D) super- 
conducting slab by a 2D simulation cell assuming the 
vortex lines are parallel to the cell edges^, in the pres- 
ence of a regular square array of pinning sites with period 
a and randomly distributed pins of comparable concen- 
tration. Pinning potential, produced by a single site, is 
modeled by a parabolic function, with U reg (U ran ) and 
o~reg (f r ati) being the depth and size of potential wells of 
regular (random) origin; a reg , a ran "C a and a ran ~ a reg . 
We will focus mostly on the weak-disorder regime, i.e., 
when the maximum pinning force by one regular site, 
f reg = 2U reg /o~ reg , is significantly larger than that for 
a random site, f ran = 2U ran /a ran . Vortices are treated 



within the London approximation (i.e., as point-like par- 
ticles). We consider the case of the first matching field, 
i.e., when the number of vortices is equal to the num- 
ber of regular pins. To study the motion of vortices, we 
use molecular-dynamics simulations, and we numerically 
integrate the overdamped equations of motior^- - 



(1) 



Here, fi is the total force acting on vortex i; f™ and f^ p 
are the forces due to the vortex-vortex and vortex-pin 
interactions, respectively; i d is an external driving force 
(i.e., a Lorentz force created by an applied current); ff is 
the thermal stochastic force. A simulation cell contains 
20 x 20 regular pins, and we use periodic boundary condi- 
tions to simulate an infinite array. To explain the results 
of numerical experiments, we also use an analytical ap- 
proach: in various situations most of the vortices remain 
pinned, so that the vortex transfer occurs through col- 
lective defects, which can be described reasonably well 
by just few parameters and by using assumptions of elas- 
ticity theory. Below we will present our results for the 
following set of parameters: a = A(T), o~ reg — 0.15a, 
a ran = 0.2a, f reg = 0.6/o, where /„ = $57327r 2 A(!r) 3 , 
numbers of regular and random sites being the same. 

De-pinning of stripes (phase I). — Very weak driving 
results in no vortex motion (pinned regime). If the driv- 
ing force Fd reaches some threshold value Fj? , part of 
vortices start to move. The vortex motion is not indi- 
vidual, since vortices travel collectively in a soliton-likc 
manner, being localized within vortex rows. Moving col- 
lective structures are just depinned stripe-like defects, 
which were predicted in the static configurations^. Each 
defect contains either one extra vortex (kink) or one va- 
cancy (antikink). Kinks propagate in the direction of an 
applied force, so that they are solitary compression waves 
of vortex row, as seen from Fig. 1 which shows kink mo- 
tion. Vacancy-based antikinks flow in the opposite di- 
rection; they can be considered as decompression waves. 
The number of vortices D in a static stripe depends on 
the balance between the vortex- vortex interaction energy 
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and the energy gained due to the displacement of vortices 
from the centers of regular pins 
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where C/^p is the critical U reg , below which the inter- 
vortex interaction dominates over the regular pinning, so 
that it becomes energetically favorable for the vortex lat- 
tice to switch globally from the square configuration to 
the so-called half-pinned phas o 10 ' 11 . In this phase, vor- 
tices in odd rows remain pinned, while vortices in even 
rows are shifted with respect to them over a/2, the pe- 
riodicity of the lattice being preserved. Eq. (2) implies 
that when approaching U^ g , the stripe becomes infinitely 
long. For long stripe, i.e., D 3> 1, we found the criterion 
for the instability of the pinned stripe, in leading order 
in l/D. The threshold occurs, when the total driving for 
D vortices in the stripe becomes equal to the resistance 
force produced by the first pinned vortex in front of the 
stripe, which prevents its motion, and the last pinned 
vortex behind the stripe, 
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Eq. (3) demonstrates that the stripe behaves nearly as 
a rigid body containing D particles and this is the ex- 
planation why it is depinned at relatively low driving 
force, Ffj 1 ' <C f reg . Fig. 2 presents typical numerically 
calculated average vortex velocity, after steady flow is 
achieved, as a function of driving, where one can clearly 
see distinct dynamical regions for a weakly disordered 
regime, while for the stronger disorder case they are 
smeared out due to chaotization. Note that in the par- 
ticular initial configuration for f ran d = freg I 6 there was 
no stripe in the simulation region, which is reflected by 
the absence of any current up to the dynamical phase 
II (curve 1 in Fig. 2). The results of Eqs. (2) and 
(3) are in good agreement with our numerical simu- 
lations: depinning of stripes is predicted to occur at 
F^P w 0.28 {D w 5, as calculated from Eq. (2)), while 
in the numerical simulations this value is around 0.23 at 

frand — freg/3. 

Kinks and antikinks can collide and annihilate. An- 
other remarkable but quite rare process in the weak dis- 
order regime process is soliton sticking by bunches of ran- 
dom pins. These two processes lead to a decay of the total 
current in time. Fig. 3 shows typical time dependences 
of the current, where two examples are addressed, when 
kinks and antikinks annihilate (curve 1) and when they 
persist (curve 1'). It is obvious, however, that in infinite 
systems all the kinks and antikinks have to disappear, 
since the total number of kinks and antikinks in a single 
row is the same. Nonvanishing motion in the weak dis- 
order regime thus appears as an artifact of a finite-size 
simulation cell with periodic boundary conditions. 
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FIG. 1: (Color online) Motion of the stripe shown by three 
snapshots. Irregular blue (light gray) spots represent posi- 
tions of random pins, red (dark gray) filled circles show posi- 
tions of vortices, and regular black open circles correspond to 
periodic pins. Dashed lines are guides to the eye indicating 
positions of defects, and arrows show the direction of their 
motion. 
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FIG. 2: (Color online) The average vortex velocity as a func- 
tion of driving for different values of the random pinning force: 

fraud = freg/6 (cUI'Ve 1), freg/3 (2), frand = freg/2 (3). 

Note that static defects in the weak and intermedi- 
ate disorder regime consist not only of stripes. Stripes 
can be just parts of clusters, some of which are large 
and branched^. Naively one may expect that these clus- 
ters act as easy-channels for vortex transfer. Instead, 
we found that low driving partially heals defects. The 
reason is that clusters basically consist of two types of 
segments^: Segments of the first kind contain no vacan- 
cies or excess vortices, in contrast with segments of the 
second type, which are nothing but stripe-like defects. 
For low driving, vortices inside segments of the first kind 
move collectively to their nearest vacant pins. Such a del- 
icate healing, however, is not possible for segments of the 
second type. Therefore, stripes do persist in the sample 
and very low driving leads to fragmentation of clusters. 
A typical driving force, which heals such defects, can be 
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FIG. 3: (Color online) Time dependence of the average vortex 
velocity for different dynamical phases at fraud = 0.2. Here 
to — 4-K 2 \(T) 4 n/&o- Curves 1 and 1' correspond to phase I at 
Fd = 0.4 and 0.2, respectively. Curves 2 and 3 show phases 
II and III for Fd = 0.52 and 0.57, respectively. 

estimated by considering an infinite chain of depinncd 
vortices. An effective pinning force for vorticies from the 
chain is created by interactions with surrounding rows of 
pinned vortices 10 . By performing a summation in recip- 
rocal space, we obtain the following estimate: 
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In the relevant range of parameters, F9 iea " <C F^. 
Generation of kink-antikink pairs (phase II). — When 

driving approaches certain critical value f kink- 

antikink pairs are created spontaneously. New kinks 
move in the direction of the drive, whereas antikinks 
propagate in the opposite direction. The correspond- 
ing dynamical phase is presented in Fig. 2 by a cusp in 
curves 1 and 2. Because of the intensive generation of 
new kinks and antikinks, average current increases sig- 
nificantly compared to the first dynamical phase, as seen 
from Fig. 3. Most of the pair generation events are trig- 
gered by moving stripes in adjacent rows, which create 
an additional washboard drive due to vortex-vortex re- 
pulsion. Such a process is shown in Fig. 4 by three snap- 
shots: one of the vortices is depinned, then it creates 
an area of vortex row compression in front of itself and 
the area of decompression behind, these two areas being 
transformed into a kink and antikink. It is easy to real- 
ize that the amplitude of the additional periodic drive is 
equal to F^ heal ^ in the limit of long kinks, so the total 
effective driving in the row, next to the moving stripe, is 



Fi e//) « F d + F^ 



However, a comparison with the 



numerical results shows that this value is too low to ex- 



is only around 0.03, 

(ii) 



plain our numerical data: p( heaL ) 

while the difference between f reg and F^ 11 ' is clearly 
disorder-dependend: it is about 0.1 for fraud = /reg/3 
and 0.05 for f reg /& (besides, there were no preexisted 
kinks within a simulation region in the latter case). The 
analysis of vortex motion patterns reveals that kink- 
antikink pairs are generated not everywhere. The reason 



is that there are some weak points, where vortices are 
additionally strongly displaced in the direction of drive 
by random pins. Concentration of these weak points de- 
pends on the external drive, i.e., on the displacement ro 
of vortices inside cores, r ~ (T reg F#^ / f reg . Depinning 
then occurs, if a pinning force by a random pin is nonzero 
in this point, which already restricts the position of the 
random pin. Another condition is that the force, pro- 
duced by the random pin at the edge of the regular pin, 
together with F^ must dominate the force by the reg- 
ular site. From these two conditions, we found a value of 
the external drive, for which an appropriate position of 
a random pin exists 

jp(II) freg fran (1 0~reg / 0~ran) j^(heal) 
1 ~t" ( fran / freg) (o~reg / 0~ran) 

Starting from Fd — F^ l:t \ a probability for the cre- 
ation of kink-antikink pairs is nonzero. We found that 
above F^ 11 ^, concentration of weak points grows as ~ 

(Fd — Fj 11 ^) 3 / 2 , so that the kink-antikink pair generation 
intensifies. Moreover, some pairs start to nucleate spon- 
taneously in weak points, without the assistance of mov- 
ing defects (in particular, for the case of / re p/6, when 
no preexisted stripes were found within the simulation 
region). Note that Eq. (5) was obtained under the as- 
sumption that a depinning is due to the action of a single 
random pin, without overlap with other ones, such over- 
laps being able to produce a larger total pinning force. 
However, events of this sort are very rare and they do not 
smear F^ 11 ^ significantly. Also, we didn't take into ac- 
count the mutual repulsion of vortices, which acts against 
the vortex depinning: according to our estimates, the re- 
sistance force is too small to noticeably change Fr 1 ' . Eq. 
(5) is in a reasonably good agreement with our numerical 
results of Fig. 2 for weak disorder, f ran d = freg/3 and 
freg I '2 (in the latter case one has to drop Fj heal ' from 
Eq. (5), since there were no preexisted stripes in our 
simulation region). In the limit of very weak disorder, 



Eq. (5) reduces to F^ 1 ' » /, 



ran - F^^ , which 

clearly shows a competition between regularity and dis- 
order, the latter factor being enhances by moving kinks 
through an additional term F^ heal \ We also notice that, 
due to a continuous generation of kink-antikink pairs, 
vortex transfer within phase II no longer decays in time, 
as in phase I. Motion occurs via groups of kinks and an- 
tikinks propagating in opposite directions and creating, 
from time to time, new defects, one of which joins the 
same group, whereas its counterpartner starts to flow in 
the opposite direction. Pairs can also be created with- 
out the assistance of previously excited stripes. In steady 
flow regime, replication of defects has to be balanced with 
their annihilation, under collisions of individual kinks 
and antikinks or their groups. 

Free proliferation of kink-antikink pairs (phase III). — 
We found third dynamical phase, which shows up as the 
last distinct region in curves 1 and 2 in Fig. 2. Care- 
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FIG. 4: (Color online) The same as Fig. 1 but with generation 
of kink-antikink pair by moving stripe, shown by snapshots. 
Initially kink flows in the bottom row, then it creates a kink- 
antikink pair in the upper row. 

ful analysis of numerical data has indicated that, in this 
phase, kink-antikink pairs start to proliferate freely not 
only in weak points: preexisted kinks immediately gen- 
erate kinks and antikinks; in their turn, these defects 
create new pairs, and so on. The area of vortex motion 
rapidly extends until it covers the whole simulation re- 
gion. Mathematically, in this regime, washboard drive 
produced by moving kink, is already strong enough to 
generate kink-antikink pairs by itself. The correspond- 
ing driving is given by 

t d ~ Jreg ~ * d ■ (O) 

In phase III, no "islands" are found, where vortices can 
be pinned for a long time. According to Eq. (6), F^ 111 ^ 
is 0.57, and this value again is in good agreement with 
the numerical results, shown in Fig. 2. 

Let us now briefly discuss the dynamics of vortices, 
when the disorder is no longer very weak, so that it not 
only triggers motion of defects, but also significantly af- 
fects it. The general tendency is that disorder smears 
out well-separated dynamical phases, as seen from Fig. 
2. Stripes can now easily jump from row to row and 
bend. However, very surprisingly, soliton-like origin of 
the vortex transfer is extremely robust, up to the regime 
of strong disorder, f ran ~ jreg- In the regime of very 
strong disorder, f ran 3> f reg , vortex flow is localized in 
narrow streams, in which vortices flow one by one. An- 
other mechanism of vortex transfer is pumping, when 
vortices are pushed into some traps, until their mutual 



repulsion breaks the blockade. 

Note that here we analysed the first-matching-field 
regime, while a little imbalance between the numbers of 
regular pins and vortices could serve as an additional 
source of disorder. As was shown in Refsj^, this im- 
balance results in different dynamical regimes includ- 
ing vortex flow in "incommensurate rows" and negative- 
differential-resistivity (NDR) parts of the Vl-curve of N- 
— ^ and S-type^. Very recently, the first experimental ob- 
servation of the N-type NDR phase has been reported 12 . 
(See also a related experiment^ on a triangular array of 
pins, where channeling of vortices can be suppressed by 
the random removal of pinning sites 14 .) 

Conclusions. — We studied the competitive effect of 
periodic square and weak random pinning potentials on 
the dynamics of vortices in two dimensions. We found 
new dynamical phases, which are established through a 
cascade of transitions. There are three phases, in which 
vortices move in a soliton-like collective structure trav- 
elling within individual vortex rows. These are kinks, 
each containing an excess vortex and moving in the di- 
rection of an external drive, and antikinks, flowing in the 
opposite direction and containing a vacancy. When col- 
liding, kinks and antikinks annihilate. Initial motion is 
triggered by depinning of preexisted static kinks and an- 
tikinks. In the second regime, moving defects generate 
secondary kink-antikink pairs in adjacent rows at certain 
weak points, which are more corrupted by disorder. In 
the third regime, these pairs are freely generated by mov- 
ing kinks and antikinks, i.e., without help of disorder. 

Although we have concentrated on vortices in super- 
conductors, it is clear that similar dynamical regimes 
will be realized in other two-dimensional systems with 
periodic lattice potentials, containing repelling particles. 
Moreover, disorder- induced kink-antikink generation un- 
der an external drive can appear as a rather universal 
phenomenon, which exists for systems and lattice poten- 
tials of various dimensionalities. 
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